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NUMERICAL  CALCULATIONS  OF  MERTEN'S  CORRECTION  FOR 
ZERO  AND  FIRST  ORDERS 

by 

Bill  D.  Cook 


The  simplified  theory  of  Raman  and  Nath  for  the  diffraction 
of  light  by  ultrasonic  waves  is  valid  only  for  a small  range  of  experi- 

i 

mental  values.  Mertens-  has  supposedly  extended  this  range  by  a per 
turbatioa  method.  It  is  the  purpose  of  this  memorandum  to  give  certain 
numerical  values  which  will  allow  the  computation  of  the  light  inten- 
sities in  the  zero  and  first  diffracted  orders  from  Merten's  results. 

The  results  of  Mertens  for  the  light  amplitude  of  the  nth 
order  can  be  expressed  as 

♦n  = Jn(v)  + i H An  (v)  + **  Bn  (v)  (l) 

where  J (v)  is  the  nth  order  Bessel  function  of  argument  v.  The 
**  2\  7T  L 

parameter  H is  where  \ is  the  wavelength  of  light,  L is  the 


path  length  of  the  light  through  the  ultrasonic  beam,  n is  the  index 

of  refraction  of  the  undisturbed  medium,  and  \ is  the  wavelength  of 
2 

sound  . An(v)  and  Bn  (v)  are  power  series  expressed  as  follows 


A = 
n 

°°  2m+n 

-ar  v 

m=0  n'° 

(2) 

and 

B * 

« . . 2m+n-2 

Z b V 

(3) 

n 

m=l 

where 

A 

( •- 1 Vn  r2tsH-n(2n+l)  ] 

n,m 

6 (Z2*^1)  (n.‘ )(m+n).' 

(4) 
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(-l)m  [6m  + (n+l)(10n  -7)]  [m  + l/6(2n2  + 3^-6)] 


n,m 

(60)  (22m  n)(m-l).'  (ra  + n-1).' 

(5) 

The  prime  over  the  summation  sign  in  Eq.  (3)  means  that  the 

sum  starts 

with  n = 2 when  n = 0.  The  light  intensity  is  given  by 

1=00*  = 
n *n  n 

(6) 

J2(v)  + H2  (A2(v)  + 2J  B (v)  ) + B2  (v)  = 

n'  n'  ' n n'  J ' n ' ' 

(7) 

/(v)  ♦ HS  Gn(v)  ♦ K4  S (v). 

where  G (v)  = A2(v)  + 2J  3 (v) 
n n n n ' 

(8) 

Mertens  gives  only  the  first  two  terms  of  Eq.  (8)  in  his  final  result. 
Tables  I and  II  give  numerical  values  needed  to  compute  1^  and  1^ 
respectively,  from  equation  (8). 

The  parameter  H is  completely  determined  by  the  experi- 
mental condition.  For  water  and  the  Hg  green  line,  it  can  be  expres- 
sed as 

H = 1.5  x 10'2  F2L  (9) 

where  F is  the  frequency  in  megacycles  and  L is  the  width  of  the  sound 
beam  in  cm. 

Table  III  gives  the  values  of  the  a and  b in  order  that 

n,  m n,  m 

one  may  calculate  A and  B for  any  value  of  v. 
n n 
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TABLE  I TABLE  II 


I 

II 

I 

II 

V 

G0(v) 

B2 

0 

G^v) 

B? 

0.2 

.000411 

_ 

-.00201 

.000017 

0.4 

.001574 

- 

-.000717 

.000061 

0.6 

.003295 

.000003 

-.001306 

.000114 

0.8 

.005894 

.000008 

-.001262 

.000156 

1.0 

.007260 

.000015 

-.001311 

.000168 

1.2 

.OO89I8 

.000023 

-.000053 

.001415 

1.4 

.010078 

.000029 

.002318 

.000086 

1.6 

.010680 

.000029 

.005777 

.000027 

1.8 

.010807 

.000022 

.010095 

- 

2.0 

.010669 

.000010 

.014855 

.000043 

2.2 

.010653 

- 

.019520 

.000184 

2.4 

.OIO697 

.000006 

.023517 

.000431 

2.6 

.011691 

.000044 

.026343 

.000763 

2.8 

.013378 

.000131 

.027660 

.001133 

3.0 

.015895 

.000280 

.027366 

.001472 

3.2 

.019095 

.000494 

.025630 

.001705 

3.4 

.022680 

.000761 

.022877 

.OOI77I 

3-6 

.026245 

.001053 

.019730 

.OOI636 

3.8 

.029348 

.001327 

.016909 

.001315 

4.0 

.031593 

.001534 

.015113 

.000872 

4.2 

.032707 

.001624 

.014898 

.000415 

4.4 

.032603 

.001564 

•016573 

.000083 

4.6 

.031410 

.001347 

.020137 

.000014 

4.8 

.029465 

.001000 

.025287 

.000316 

5.0 

.027269 

.000592 

.031354 

.001038 

5.2 

.025404 

.000222 

.037605 

.002142 

5.4 

.024443 

.000013 

.043166 

.003511 

5.6 

.024846 

.000087 

.047262 

.004947 

5-8 

.026880 

.000545 

.049339 

.006205 

6.0 

.030572 

.001440 

.049161 

.007040 

6.2 

.035253 

.002650 

.047068 

.007326 

6.4 

.040874 

.004187 

.043270 

.006879 

6.6 

.046559 

.005826 

. 038682 

.005775 

6.8 

.051561 

.007332 

.034261 

.004198 

7.0 

.055202 

.008458 

.030973 

.002459 

7.2 

.057001 

.008965 

.029614 

.000952 

7-4 

.056756 

.008772 

.030706 

.000086 

7-6 

.054604 

.007798 

.034356 

.000205 

7-8 

.051008 

.006184 

.040250 

.001513 

8.0 

.046697 

.004193 

.047667 

.00  to  15 
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TABLE  I (continued) 


V 

G0(v) 

B2 

0 

8.2 

.042558 

.002205 

8.4 

.039492 

.000664 

8.6 

.038265 

.000004 

8.8 

.039372 

.000568 

9.0 

.042943 

•002533 

9.2 

.048707 

.005853 

9,4 

.056014 

.010237 

9 6 

.063935 

.015169 

9.8 

.071399 

.019971 

10.0 

.077363 

.023906 
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Demonstration  of  the  “Least  Stable  Waveform” 
of  Finite  Amplitude  Waves 

M.  A.  Breazeale  and  W.  W.  Lester 
Depart  mini  of  Physics,  Michigan  Stoic  University , 

East  Lansing,  Michigan 
(Received  July  28,  1961) 

A backward-sloping  ultrasonic  wave  is  generated  by  reflection 
of  a finite  amplitude  distorted  wave  from  a pressure-release  bound- 
ary. The  difference  between  the  behavior  of  this  wave  and  that  of  a 
wave  distorted  in  the  usual  sense  is  demonstrated. 

AS  is  well  known,  an  ultrasonic  wave  of  large  amplitude  can 
undergo  distortion  so  that  the  slope  of  its  leading  edge  is 
much  greater  than  that  of  its  trailing  edge.  If  such  a nonsymmetric 
wave  is  reflected  from  a boundary,  two  extreme  cases  are  possible: 
The  waveform  may  be  unchanged,  as  in  the  case  of  a perfectly 
rigid  reflector,  or  it  may  be  inverted,  as  in  the  case  of  a perfect 
pressure- release  reflector.  For  a pressure-release  reflector,  the  re- 
sulting waveform  would  have  the  slope  of  its  trailing  edge  greater 
than  the  slope  of  its  leading  edge.  Such  a waveform  has  been 
described  by  Fay1  as  the  “least  stable  waveform”  in  contrast 
with  the  more  usual  case  of  the  most  stable  waveform.  In  the 
least  stable  waveform  the  fundamental  harmonic  component  can 


Fig.  1.  Effect  of  phase 
shift  of  second  harmonic  in 
resonant  transducer,  (a)  In- 
cident distorted  wave,  (b) 
resulting  CRO  wave. 


4—  (a) 


(b) 


increase  with  distance,  while  the  higher  harmonics  decrease  with 
distance.2  Thus,  an  initially  distorted  wave  can  become  undistorted 
as  it  progresses.  The  following  experiment  was  designed  to  demon- 


strate the  inversion  of  waveform  on  reflection  from  a pressure- 
release  boundary  and  this  decrease  of  distortion  with  distance. 

A pi ’ssure-release  reflector  was  constructed  by  stretching  a 
plastic  membrane  over  an  air  chamber,  and  a rigid  one  by  mount- 
ing an  aluminum  plate  so  it  could  be  easily  attached  in  place  of 
the  pressure  release  boundary. 

The  experiment  was  performed  with  high  intensity  2-Mc  ultra- 
sonic pulses  in  water.  The  waveform  of  the  ultrasonic  pulse  at 
various  distances  before  and  after  it  had  been  reflected  was  moni- 
tored by  use  of  a barium  titanate  transducer.  Since  the  receiver 
resonated  near  the  second  harmonic,  this  harmonic  was  accentu- 
ated in  the  waveform  displayed  on  an  oscilloscope.  Further,  the 
receiver  resonance  caused  a phase  shift  of  almost  90°  in  the  second 
harmonic  component  displayed  on  the  oscilloscope.  This  made  it 
possible  to  detect  very  small  amounts  of  second  harmonic  by  ob- 
serving the  "flattening”  of  the  bottom  of  the  oscilloscope  trace 
and  the  “sharpening”  of  the  top.  This  is  illustrated  in  Fig.  1.  For 
a distorted  wave  such  as  that  illustrated  in  (a)  we  get  the  oscillo- 
scope trace  shown  in  (b)  when  the  phase  shift  of  the  second  har- 
monic is  90°.  The  increase  of  waveform  distortion  with  distance 
could  be  observed  by  watching  this  type  of  distortion  on  the  scope. 
Since  there  is  a phase  shift  of  180°  in  each  of  the  harmonics  when  a 
rigid  reflector  is  used  and  0°  when  a presssure-release  reflector  is 
used,  the  corresponding  oscilloscope  traces  differed  markedly  from 
each  other. 

Figure  2 gives  a series  of  pictures  of  oscilloscope  traces  which 
show  the  transducer  output  at  increasing  distances  from  the  source. 
It  can  be  seen  that  the  distortion  of  the  waveform  increases  with 
distance.  As  indicated,  reflection  from  a solid  boundary  produced 
the  waveforms  shown  in  the  second  row  of  pictures  where  it  can 
be  seen  that  the  distortion  continues  to  increase  with  increasing 
distance.  On  the  other  hand,  the  lower  part  of  the  figure  show's  that 
the  wave  on  reflection  from  a pressure-release  boundary  is  inverted ; 
i.e.,  it  is  distorted  in  the  “wrong  direction.”  Therefore,  the  distor- 
tion decreases  with  increasing  distance  until  at  a distance  from 
the  pressure-release  reflector  almost  equal  to  that  between  the 
transducer  and  the  boundary,  the  waveform  is  again  essentially 
sinusoidal.  On  progressing  farther,  the  wave  becomes  sinusoidal 
and  then  distorts  again.  Considerations  of  this  type  of  distortion 
might  be  necessary  if  one  were  concerned  with  reflection  of  finite 
amplitude  waves  from  the  free  surface  of  liquids. 

The  authors  wish  to  express  their  appreciation  to  Professor 
E.  A.  Hiedemann  for  the  interest  he  has  shown  in  this  work,  and 
to  the  Office  of  Naval  Research  for  their  sponsorship. 

> R.  D.  Fay,  J.  Acoust.  Soc.  Am.  29,  1200  (1957). 

1 Isadore  Rudnick,  J.  Acoust.  Soc.  Am.  30.  564  (1958). 
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Fig.  2.  Distortion  of  a finite  amplitude  wave  as  it  progresses  from  a sinusoidally  vibrating  souice  and  is  reflected  from  a rigid 
boundary  or  a pressure  release  boundary.  The  boundaries  are  90  cm  from  the  source. 
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An  optical  method  is  used  to  measure  the  energy  ratio  of  reflected  and  incident  ultrasonic  waves  at  a 
liquid-solid  interface.  The  ultrasonic  velocities  in  the  solid  are  calculated  from  the  angles  of  maximum  re- 
flection in  the  liquid. 


THE  intensity  ratio  of  reflected  to  incident  ultra- 
sonic waves  at  a liquid-solid  boundary  as  a func- 
tion of  angle  of  incidence  is  given  by  Ergin1  as 


(tH 


cos/3— A cosa(l  — B)"|2 
cos/3 + A cosa(l  — B)\ 


where  a is  the  angle  of  incidence  in  the  liquid  measured 
from  a line  normal  to  the  interface,  /3  and  y are  the 
angles  of  refraction  of  the  longitudinal  and  shear  wave 
in  the  solid.  The  quantities  A and  B are  defined  by 

A = piV  l/piVi,  (2) 

B=  2 siny  sin2y[cosy—  (Fa/ Ft)  cos/8],  (3) 

where  pi  and  p2  are  the  densities  of  the  liquid  and  the 
solid,  respectively,  and  V t is  the  velocity  of  the  incident 
wave  in  the  liquid ; V l and  Fa  are  the  velocities  of  the 
refracted  longitudinal  and  shear  waves  in  the  solid. 

Substituting  accepted  values  for  the  densities  and 
velocities  of  water  and  Plexiglas  in  Eq.  (1)  one  obtains 
the  curve  shown  in  Fig.  1(a)  for  the  intensity  ratio 
(. R/I)2 . Figure  1(b)  shows  this  ratio  for  a water- 
aluminum  boundary.  The  ultrasonic  wave  is  incident 
in  the  water  in  both  cases.  In  order  to  obtain  these 
curves  one  has  to  use  the  appropriate  angles  /3  and  y 
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for  a given  angle  of  incidence  a.  These  angles  are  found 
from  Snell’s  law 

F//Ft=sina/sin/3,  F//Fs  = sina/siny.  (4) 

It  can  be  seen  from  Eq.  (4)  that  sina=  Fj/F t or 
sina=  Vi/V s at  the  critical  angles  for  the  longitudinal 
and  shear  wave  where  sin/3  or  siny  equal  unity.  One 
can  obtain  F i and  F s if  Vj  is  known  provided  the 
critical  angles  can  be  located.  Equation  (1)  and  Fig.  1 
show  that  the  ratio  (R/I)2— l at  the  critical  angles. 
The  associated  peaks  in  the  intensity  of  the  reflected 
wave  can  be  seated  experimentally  and  can  be  used 
to  calculate  Vl  and  Fa  for  the  solid.1 

An  optical  method  is  used  to  find  the  angle  of  in- 
cidence at  which  the  intensities  of  the  reflected  and 
incident  beams  are  equal.  The  arrangement  is  shown 
in  Fig.  2.  The  solid  sample  and  the  transducer  are 


PhoUmultipifer 


Fio.  1.  Intensity  ratio  of  reflected  to  incident  wave  as  a func- 
tion of  angle  of  incidence  for  (a)  a water- Plexiglas  boundary 
where  Vs<Vi<Vl  and  (b)  a water-aluminum  boundary  where 
V,<V.,<Vl. 

* Presented  at  the  62nd  meeting  of  the  Acoustical  Society  of 
America,  Cincinnati,  November  1961. 

t NSF  Undergraduate  Research  Paaicipant. 

' K.  Ergin,  Bull.  Seism.  Soc.  Am.  42,  349  (1952). 


Fig.  2.  Diagram  of  the  optical  arrangement. 

placed  in  a tank  filled  with  water.  While  the  angle  of 
incidence  is  changed  by  rotating  the  transducer,  the 
sample  is  also  rotated  in  such  a manner  that  the  re- 
flected sound  beam  remains  at  right  angles  to  the 
collimated  light  beam.  The  reflected  ultrasonic  wave 
produces  a diffraction  pattern  in  the  plane  of  the  pho- 
tomultiplier. The  light  intensity  in  the  nth  order  of  the 
diffraction  pattern  is  given  by 

/»=•/,’(«-),  (5) 

where  v is  proportional  to  the  amplitude  of  the  ultra- 
sonic wave  producing  the  diffraction  pattern.  Keeping 
the  output  of  the  transducer  constant  and  measuring 
the  light  intensity  in  the  zero  order,  one  finds  a pro- 
nounced dip  at  that  angle  of  incidence  where  the  re- 
flected ultrasonic  wave  is  most  intense.  Figure  3(a) 

* W.  G.  Mayer,  J.  Acoust.  Soc.  Am.  32,  1213  (1960). 
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shows  the  zero  order  light  intensity  for  an  8-Mc  con- 
tinuous ultrasonic  wave  reflected  from  a water-Plexi- 
glas  boundary.  For  convenience  the  intensity  of  the 
incident  wave  is  kept  low  enough  so  that  the  zero-order 
Bessel  function  is  positive  for  all  possible  values  of  v 
of  the  reflected  wave.  The  critical  angle  is  located  at 
33.5°  from  which  one  finds  Ft=2700  m/sec,  using 
Eq.  (4).  The  velocity  in  water  Vj  can  be  determined 
by  measuiing  the  spacing  between  the  lines  of  the 
diffraction  pattern  produced  by  the  reflected  wave  in 
the  liquid.  Figure  3(b)  shows  the  intensity  of  the  re- 
flected wave  obtained  from  the  data  given  in  Fig.  3(a). 
The  theoretical  curve  predicted  by  Eq.  (1)  is  also  given. 

The  same  technique  is  used  to  measure  Us  if  F.s>  l'/. 
In  this  case  the  light  intensity  in  the  zero  order  reaches 
a minimum  at  the  critical  angle  for  the  shear  wave  and 
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Fig.  3.  Critical  angle  for  longitudinal  wave  in  Plexiglas:  (a) 
zero-order  light  intensity  in  diffraction  pattern  produced  by  re- 
flected wave  in  water;  (b)  corresponding  values  of  (/?//'’.  Solid 
line  shows  theoretical  values. 
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ANGLE  OF  INCIDENCE 

Fig.  4.  Critical  angle  for  shear  wave  in  glass:  (a)  zero-order 
light  intensity  in  diffraction  pattern  produced  by  reflected  wave 
in  water;  (b)  corresponding  values  of  (/?//)*.  Solid  line  shows 
theoretical  values. 

remains  at  that  level.  An  example  is  given  in  Fig.  4(a) 
which  shows  the  measured  light  intensity  in  the  zero 
order  produced  by  a wave  reflected  from  a water-glass 
boundary  at  angles  in  the  vicinity  of  the  critical  angle 
for  the  shear  wave.  The  corresjKmding  intensity  of  the 
reflected  wave  is  shown  in  Fig.  4(b).  The  critical  angle 
for  the  shear  wave  is  25.4°  corresponding  to  a shear 
wave  velocity  of  3445  m sec. 

It  should  be  noted  that  this  analysis  does  not  include 
surface  waves  or  plate  transmission  phenomena.  The 
method  given  here  has  the  advantage  that  the  velocity 
of  the  longitudinal  and  shear  wave  in  the  solid  can  be 
calculated  without  having  to  observe  the  waves  in  the 
solid  directly. 
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A theory  is  developed  for  the  diffraction  of  light  by  two  spatially  separated 
parallel  ultrasonic  progressive  waves  of  different  frequency.  The  preliminary  theo- 
ries of  Raman  and  Nath  [C.  V.  Raman  and  N.  S.  Nath,  1’roc.  Indian  Acad.  Sci. 
A 2,  406—412;  413  — 420  (1935)]  for  normal  and  oblique  incidence  are  taken  to  be 
valid.  The  resulting  equations  are  extensions  of  earlier  results  of  R.  Mertens, 
Z.  Physik  160,  291—296  (I960).  The  predicted  periodicity  of  the  diffraction  spec- 
trum with  increasing  sound  beam  separation  agrees  with  the  well  known  periodicity 
of  the  light  intensity  distributions  in  the  Fresnel  zone  of  the  phase  grating  formed 
by  the  first  ultrasonic  wave.  Results  of  numerical  calculations  are  presented  to 
illustrate  features  of  the  theoretical  results,  as  reflected  in  the  first  order  of  diffrac- 
tion for  3-0  and  6-0  Me  ultrasonic  waves  in  water. 


Introduction 

Theoretical  results  have  been  given  by  Rao1,  Murty*,  and  Mertens3 
for  diffraction  of  a wide  light  beam  by  an  ultrasonic  wave  consisting 
of  two  commensurable  frequency  components  with  arbitrary  relative 
phase.  These  theoretical  results  are  for  simultaneous  diffraction  of  light 
by  the  two  frequency  components  contained  in  the  same  sound  beam. 
Mertens4  recently  pointed  out  that  simultaneous  diffraction  and 
successive  diffraction  by  two  separate  parallel  ultrasonic  beams  are  not 
the  same.  Murty  and  Rao5  obtained  very  good  agreement  between 
their  experimental  results  from  a successive  diffraction  expo  "iment  and 

1 Rao,  B.  R.:  I’roc.  Indian  Acad.  Sci.  A 29,  16—27  (1949). 

3 Murty,  J.  S. : J.  Acoust.  Soc.  Amer.  26,  970  — 974  (1954). 

3 Mertens,  R.:  Proc.  Indian  Acad.  Sci.  A 8,  288  — 306  (1958). 

4 Mertens,  R.:  Z.  Physik  160,  291—296  (I960). 

5 Murty,  J.  S.,  and  B.  R.  Rao:  Z.  Physik  157,  189  — 197  (1959).  The  agreement 
obtained  is  indeed  surprising,  as  it  is  doubtful  that  progressive  waves  were  obtained 
at  3.0  and  1.5  Me  in  a metal  tank  only  12  in.  long. 
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values  calculated  from  simultaneous  diffraction  theory.  Also,  the  dif- 
fraction effects  of  two  parallel  ultrasonic  waves  have  proven  useful,  in 
a limited  range,  in  investigations 8,7  of  finite-amplitude  distortion:  Simul- 
taneous diffraction  theory  was  used.  But,  as  Mertens  has  stated, 
simultaneous  and  successive  diffraction  differ  negligibly  for  the  cases 
considered  experimentally  by  Murty  and  Rao.  In  a limited  range, 
simultaneous  and  successive  diffraction  spectra  are  indistinguishable. 
Hargrove8  made  some  preliminary  experimental  studies  of  successive 
diffraction  and  reported  qualitative,  but  lack  of  quantitative,  agreement 
between  his  experimental  results  and  simultaneous  diffraction  theory. 
Mertens’  4 paper,  giving  expressions  for  the  amplitudes  of  light  diffracted 
by  two  parallel,  adjacent  (no  space  between),  ultrasonic  beams  with 
integer  frequency  ratio  and  arbitrary  relative  phase,  appeared  shortly 
after  the  experimental  results  of  Hargrove  were  reported.  The  present 
theoretical  development  was  undertaken  to  include  the  spatial  separation 
of  ultrasonic  beams  which  existed  in  the  experimental  arrangement  of 
Hargrove. 

Raman  and  Nath8’10  developed  a theory  to  explain  the  diffraction 
of  light  by  sound  for  the  case  of  a wide  plane  wave  of  light  passing  through 
a plane  sinusoidal  ultrasonic  wave  at  normal  incidence.  The  theory 
predicts  that  the  light  is  diffracted  at  discrete  angles  given  by 

sin  = — rA/A*,  (1) 

where  r is  zero  or  a positive  or  negative  integer,  and  A and  A*  are  the 
wavelengths  of  light  and  sound,  respectively.  Equation  (1)  is  valid  for 
any  plane  periodic  sound  wave.  The  Raman-Nath  theory  predicts  that 
for  a progressive  sinusoidal  ultrasonic  wave  the  normalized  intensity  I, 
in  the  r-th  order  of  diffraction  is 

/,=/,».  (2) 

where  the  Raman-Nath  parameter  v is  approximately  proportional  to 
the  sound  pressure  amplitude  and  given  by 

v=2n  fiL/A.,  (3) 

where  ji  is  the  peak  change  in  refractive  index  caused  by  the  sound 
pressure  and  L is  the  width  of  the  sound  beam.  Equation  (2)  is  valid 
if  the  light  wavefront  can  be  considered  to  be  changed  only  in  relative 
phase  as  it  passes  through  the  sound  beam.  This  assumption  has  been 

9 Zankel,  K.  L.:  J.  Acoust.  Soc.  Amer.  32,  707—713  (I960). 

7 Mayer,  W.  G.,  and  E.  A.  Hiedemann:  J.  Acoust.  Soc.  Amer.  32,  706— 70S 
(I960). 

9 Hargrove,  L.  E. : J.  Acoust.  Soc.  Amer.  32,  940  (A)  (I960). 

• Raman,  C.  V.,  and  N.  S.  Nath:  Proc.  Indian  Acad.  Sci.  A 2,  406—412  (1935) • 

10  Raman,  C.  V.,  and  N.  S.  Nath:  Proc.  Indian  Acad.  Sci.  A 3,  75  — 84  (1936). 
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considered  justifiable11- 12  for  conditions  under  which 

(2  nLtylfaW&N,  (4) 

where  fi9  is  the  refractive  index  of  the  undisturbed  medium  and  N lies 
in  a range  1 ^ IV  ^ 2,  depending  on  the  accuracy  required. 

Proceeding  from  assumed  validity  of  the  Raman-Nath  theory,  we 
shall  obtain  expressions  for  the  diffraction  of  light  by  two  spatially 
separated  parallel  ultrasonic  waves  of  different  frequency.  This  devel- 
opment follows  the  method  used  by  Mertens4  and  extends  his  results 
to  more  general  cases.  Some  results  of  numerical  calculations  will  be 
presented  to  illustrate  the  dependence  of  the  first  diffraction  order  light 
intensity  on  various  parameters. 

Development  of  the  Theory 

Consider  two  parallel  ultrasonic  beams  with  width  Lm  and  L„  sepa- 
rated by  a distance  V as  indicated  in  Fig.  1 . Let  plane  monochromatic 


Fig.  1.  Coordinate  axis  for  the  case  of  two  spatially  separated  parallel  ultrasonic  beams  and  the  schematic 
path  of  a typical  light  ray 

light  be  incident  in  the  + 2-direction.  The  first  sound  beam  produces 
a change  in  refractive  index  given  by 

l*m  (*.  =i«o  + t>m  s*n  2?r  (i m v*t  — m x/A*  + dj2jt)  (5) 

where  r*  is  the  ultrasonic  fundamental  frequency  and  dm  is  the  relative 
phase  of  the  w-th  harmonic.  Similarly,  for  the  second  sound  beam, 

=/i0+fi„sin2x(nv*t  — nxll*+dj2Ji).  (6) 

11  Extermann,  R.,  and  G.  Wannier:  Helv.  phys.  Acta  9,  520—532  (1936). 
**  Rytov,  S.  M. : Diffraction  de  la  lumiire  par  les  ultra-sons.  Actuality  scienti- 
fiques  et  industrielles,  613.  Paris:  Hermann  & Cie.  1938. 
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We  have  assumed  that  there  exists  a fundamental  frequency  v*  such  that 
« and  m are  integers.  The  enumeration  of  diffraction  orders  and  the 
corresponding  diffraction  angles  shall  be  those  pertaining  to  the  funda- 
mental frequency. 

According  to  the  preliminary  theory  of  Raman  and  Nath,  the  light 
amplitude  at  the  point  (*',  Ln)  is  given  by 

?(*'.  TJ=exp[~2";^]X 

X exp  [—  tvmsin  2n(mv*t  — fnx'/k*  -f  dj 2n\  exp  (2 Jr  i vt) , 

where  v is  the  frequency  of  the  incident  light.  From  the  identity 

+ 00 

exp  (—  »<jsin  b)  = exp  [ia  sin  (— b)]  = £ Jp(a)  exP  {—ipb) , (8) 

P——oo 

Eq.  (7)  becomes 

r . , i +oo 

<p(x',  Lm)  = exp | exp (2n ivt)  £ Jp(v Jx 
X exp  [2  nip  (mx'/A*  — mv*t  — bj2n )] . 


Consider  now  the  p-th  term  of  Eq.  (9)  which  represents  light  propagating 
in  the  direction  &mp  = — sin-1  [mp  klfi0A*]  with  amplitude  Jp(vm). 
This  light  component  progresses  to  the  plane  z = Lm+  L'  through 
medium  undisturbed  by  ultrasonic  waves.  However,  the  light  path 
(see  Fig.  1)  is  displaced  from  x'  to  x" . At  the  point  (x",  Lm  + L') 
the  amplitude  of  the  />-th  component  is 


*<*"• L ■ + V)  - «p  pHss?-]  “P  [■  x 

X exp  [2 Jt  i (v  ~mpv*)  <]  exp  [ exp  (-  ip  6 m)  Jp(vm) . 


(10) 


Equation  (10)  expresses  the  amplitude  of  the  plane  wave  of  light  incident 
on  the  second  sound  beam,  making  an  angle  &mp  with  the  2-axis.  Using 
the  Raman-Nath  elementary  theory  for  oblique  incidence13,  Eq.  (10) 
must  be  multiplied  by 


exp 


Lncosf)mp 

(— 2n.il k)  / ftn(s,t)ds 

0 


(11) 


where  the  integral  in  (11)  represents  the  optical  path  length  of  the  p-th 
component  in  the  second  sound  beam.  The  coefficient  of  this  integral 
should  be  — 27ii  {v  — mpv*))c  but  is  here  approximated  by  —2n ijk 


13  Raman,  C.  V.,  and  N.  S.  Nath:  Proc.  Indian  Acad.  Sci.  A 2.  413  — 420  (1935). 
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since  »'<»  and  only  a few  terms  of  Eq.  (9)  are  significant.  From 
Eq.  (2)  we  obtain 

«„  (s.  <)  =i“o  + sin  2 n [n  v*t-n(x-s  sin  &„p)fX*  + Sj2n] . (12) 
The  intregal  in  (11)  then  becomes 


L*/COS  p 

f f<n(s,t)ds=  ^ 


yi  7i sin  vmp 


X sin  2 nn(v*t-  ^ tan  dm pJ  + <5„  sin  " tan  dm t 


The  value  of  the  integral  is  substituted  into  (11)  and  this  result  mul- 


tiplies Eq.  (10).  The  result  of  this  multiplication  is  expanded  (q  = sum- 
mation index)  according  to  Eq.  (8).  The  17-th  component  represents  the 
amplitude  of  a final  light  component  which  has  been  successively  dif- 
fracted d„p  by  the  first  sound  beam  and  0„q  by  the  second  sound  beam, 
and  is  given  by 


9V,~/»<0/t{— [^tan*.,]}  X 

xe*p[— -t-?aj]exp[^H^-(Z.„  + Z.')tan0H,|x  . (14) 
X exp[  =igia. 


Time  and  space  terms  which  do  not  contribute  to  the  final  intensities 
have  been  deleted  from  Eq.  (14).  In  the  argument  of  the  last  exponen- 
tial factor  (with  L'=  0),  Mertens4  essentially  approximated  1/cos  &mp 
by  unity  and  this  term  no  longer  contributed  to  the  final  intensities. 
However,  as  this  last  exponential  factor  accounts  for  the  different  optical 
phases  of  the  light  components  which  combine  to  form  a given  order 
of  the  final  spectrum,  we  must  retain  it.  The  different  phases  arise  from 
the  different  (for  different  p ) path  lengths  from  the  plane  z=Lm  to  the 
plane  z — Lm  + Ln  + L'  in  a medium  with  refractive  index  //„ . The  next 
to  last  exponential  factor  in  Eq.  (14)  expresses  the  effect  of  the  average 
ultrasonic  phase  difference  (i.e.,  in  addition  to  dm  and  6„)  along  the 
oblique  light  path  through  the  second  sound  beam.  The  second  exponen- 
tial factor  expresses  the  effect  of  the  additional  ultrasonic  phase  differ- 
ence between  the  points  (x',Lm)  and  (x,  Lm  + Ln  + L')  of  emergence 
of  light  from  the  two  ultrasonic  beams.  The  effect  of  the  actual  phase 
differences  of  the  two  ultrasonic  waves  is  expressed  by  the  first  exponen- 
tial factor  in  a manner  identical  with  that  for  simultaneous  diffraction. 
The  factors  multiplying  vn  in  the  argument  of  the  Bessel  function  of 
order  q give  an  effective  Raman-Nath  parameter  for  oblique  incidence. 

The  r-th  order  of  the  final  spectrum  with  frequency  (v  — r v”)  makes 
an  angle  &r  =&mp  +dnq  with  the  2-axis.  Summing  all  <pPq  such  that 
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mp  -\-nq=r  we  obtain  the  amplitude  in  the  r-th  order  to  be 

<Pr=  2'  <PP « 05) 

M“~  00 

where  the  prime  on  the  summation  indicates  summation  on  all  integer 
values  of  p and  q such  that  mp  + nq  =r,  r=  0 or  positive  or  negative 
integers. 

For  the  special  cases  m = 1 and  dm  = 0 or  n = 1 and  d„  = 0,  the  sum- 
mations can  be  written  in  a straightforward  manner  with  a single 
summation  index  and  no  prime  on  the  sum. 

For  m = 1,  p=r  — nq  (always  an  integer)  giving 

M Sin  [lT  ‘antf,_„,j}  X 

X exp(-t>d„)exp|— — (Ln  + L')  tan  f?,_„,Jx  05') 

x exp  [^(I./2)  exp  [ ^ . 

For  n = \,  q =r  — mp  (always  an  integer)  giving 

*' = Xj” (Vn)  J,mp  sin  tan  *4x 

X exp(  ip  <y  expl + L')  tan  &m/>\  x 

L * J (15  ) 

X exp  j ~ 2 * m p)  (/„„/2)  tan  &mp\  x 

X exp  [ — ^ i~—  ---**  t-7-1 . 

[ A coa»mp  I 

Equations  (15')  and  (15")  with  L'=  0 correspond  to  Mertens4  Eqs.  (9) 
and  (11)  respectively. 

We  may  express  <pr  in  approximate  forms  which  are  more  suitable 
for  numerical  calculations  by  using 

dk  sb  tan  dk  sin  ■&k  — — k A/(//0  A*) , (16) 

and  the  variables  1 /cos  0k  = sec  ■&k  sb  1 + #*/2,  (17) 

Q = and  V-Z&U.  (18) 

Use  of  (16),  (17),  and  (18)  gives  the  approximate  form  of  Eq.  (15)  to  be 
<Pr  = X JPK)J,  k S‘i \nJloQ  1 exp  *^Sm  + qd^  X 

P.1--* o 1 tnmpQ  j (19) 

X exp  \L%mpnqQ  + P*(Q  + Q')\. 
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For  m=\  the  approximate  form  is 

<fr ='2j,~n,  k)  j,  [»« exp  (-  *?<*..)  * 

X exp  | --  i n q (r  - nq)  Q + ± i (r  - nq)2  (Q  + Q')  j . 
For  n = 1 the  approximate  form  is 

?>r=^2  Jp{Vm)J,-mp[v i ] exP(~  X 

X exp  |-i  i m p (r-mp)  Q+^im2p2(Q  + (>')] . 


09') 


(19") 


The  more  useful  results  are  summarized  in  these  last  three  equations. 
The  light  intensities  are  obtained  from 


k I 


(20) 


Discussion 


The  inclusion  of  the  space  L'  between  the  two  sound  beams  is  a 
newly  added  parameter  in  the  theory  of  successive  diffraction.  In  most 
experimental  arrangements  there  will  be  some  beam  separation  necessary 
to  accomodate  the  crystal  mountings.  Let  us  first  consider  the  effect  on 
diffraction  of  the  beam  separation  L' . From  Eq.  (19)  we  see  that  the 
dependence  of  light  amplitude  on  L'  is  expressed  by  the  factor 


exp  (\im2pl  Q') 


(21) 


and,  being  independent  of  n and  q,  (21)  depends  only  on  L'  and  para- 
meters pertaining  to  the  first  ultrasonic  beam.  The  character  of  the 
diffraction  spectrum  should  be  periodic  as  the  argument  of  (21)  changes 
by  integral  multiples  of  2 n for  all  values  of  p.  By  expressing  the  argu- 
ment of  (21)  in  terms  of  L',  taking  p = 1,  and  equating  the  result  to 
2 Pti,  where  P is  a positive  integer,  we  obtain  the  periodicity  rela- 

tlOIlShlP  L'P  =2-P(/i0A**/m*A).  (22) 


Using  ?.*  =X*l*n  we  now  note  that 


L'P  = 2P(^}.litX)  =2PD 


(23) 


where  D is  the  distance  from  z = Lm  in  which  the  phase  modulated  wave- 
front  at  z =Lm  reappears,  as  predicted  by  Nath14  for  progressive  waves. 
At  z = Lm  + D the  light  wavefront  is  the  same  as  at  z = Lm  except  for 
being  shifted  by  AJJ/2  in  the  ^-direction.  Thus  we  see  that  the  effect 

14  Nath,  N.  S.:  Proc.  Indian  Acad.  Sci.  A 4,  262—274  (1936). 
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of  L'  equal  to  odd  multiples  of  D is  the  same  as  changing  by  i n 
radians  with  L' = 0.  Furthermore,  the  diffraction  spectra  for  any  L' 
plus  integral  multiples  of  2D  are  identical.  These  periodicity  features 
will  be  illustrated  with  numerical  examples. 

Numerical  calculations  of  the  positive  first-order  light  intensities 
were  made  using  an  electronic  computer.  Several  of  the  parameters 
were  chosen  identical  with  those  for  which  one  of  the  authors  has  experi- 
mentally investigated  simultaneous16  and  successive®  diffraction.  The 
fixed  parameters  used  for  calculations  are  v*  =)-0  Me,  A*  =0-5  mm 


(water),  A = 5461  A,  m=  2,  n=  1,  /<„  = 1-33  (water),  Lm=Ln  = 2-0 cm 
and  vt  =2-40.  The  equation  used  for  calculation,  obtained  from  Eq. 
(15").  was 

I -f-2  2 

^ exP  [2*>2<?'  -F  «?  — d2)  I . (24) 

Equation  (24)  was  evaluated  for  a2  = from  0 to  0-200  in  intervals 
of  0-025,  X = (Q  — <J2)  from  0 to  2n  in  intervals  of  .t/16,  and  Q'  from  0 to 
n in  intervals  of  jr/16.  For  a%  — 0-200  calculations  were  also  made  for 
Q'  from  0 to  ji/4  in  intervals  ji/64. 

As  one  varies  the  relative  phase  between  the  two  sound  beams,  the 
intensities  in  the  diffraction  orders  oscillate.  Our  attention  shall  be 
limited  to  amplitude  and  phase  of  the  oscillations  of  the  first-order  light 
intensity  and  the  effect  on  them  of  the  ratio  at  and  the  sound  beam 
separation  as  expressed  by  the  variable  Q'. 

Typical  first-order  light  intensity  vs  the  variable  X is  shown  in  Fig.  2. 
Note  that  the  extrema  of  the  intensity  in  this  case  do  not  occur  for  X 

16  Hargrove,  L.  E.,  and  E.  A.  Hiedemann:  J.  Acoust.  Soc.  Amer.  33,  1747—1749 
(1961). 
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Fig.  3.  The  values  of  X corresponding  to  extreme 
values  of  first-order  light  intensity  vs  (T,  for 
aj  =*  0-200.  The  open  circles  indicate  exact 
values  at  Q'  equal  to  integral  multiples  of  n/4, 
independent  of  a. 


Fig.  5.  Extrema  of  first-order  light  intensity  vs  the 
ratio  a,  for  various  values  of  Q’ : (A)  Q’  = k n/2, 
(B)  Q’  =knl2±nl\6,  (C)  Q'  = k n/2  ± n/S,  (D)  (>' 
= knl2  ± 3n/l6,  (E)  Q'=  kn/2  ± n/ 4,  where  k is  zero 
or  integral  such  that  Q'^0 


Fig.  4.  Extrema  of  first-order  light  intensity  with  respect  to  variation  in  relative  phase  and  first -order  light 
intensity  for  fixed  relative  phase,  rs  Q’,  0,  = 0-200. 
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(nor  for  d2)  equal  to  integral  multiples  of  ji,  and  that  the  curve  is  not 
symmetric  about  the  extrema.  For  the  case  of  simultaneous  diffraction, 
the  relative  phases  for  extrema  are  integral  multiples  of  n.  For  the  case 
of  successive  diffraction,  the  relative  phases  for  extrema  depend  on 
both  Q and  Q'  and  also,  in  general,  on  the  amplitudes  of  the  separate 
sound  beams.  For  the  special  cases  Q'  equal  to  integral  multiples  of 
.t/4,  the  phases  for  extrema  appear  to  be  independent  of  a2 . Figure  3 
shows  the  values  of  A'  corresponding  to  extrema  of  the  first-order 
light  intensity,  with  respect  to  variation  of  X,  vs  Q'  for  a2  =0-200. 
The  curves  for  XMax  and  XMin  for  Q'  not  equal  to  integral  multiples 
of  ji/4  are  approximate  because  the  location  of  extrema  was  estimated 
from  calculations  made  in  discrete  intervals  of  the  appropriate  variable. 
Symmetry  of  numerical  values  about  the  extreme  values  for  Q'  equal 
to  integral  multiples  of  tt/4  indicated  that  these  calculated  values 
are  truly  extrema  and  that  the  intensity  vs  X curves  are  symmetric 
about  the  extrema  for  these  special  values  of  Q'. 

Figure  4 shows  the  extrema  of  first-order  light  intensity  with  respect 
to  variation  of  relative  phase,  for  a2  =0-200,  vs  Q' . This  clearly  illus- 
trates the  periodicity  of  the  diffraction  with  beam  separation.  Another 
illustration  of  the  periodicity  is  also  shown  in  Fig.  4.  The  first-order 
light  intensity  for  constant  A anda2  =0-200  vs  Q'  is  shown  for  A = 0 and 
for  X = ±jr  (i.e.,  for  d2  — Q and  d2  = Q±7i).  Note,  in  Fig.  4,  the 
difference  between  extrema  (with  the  corresponding  changes  in  relative 
phase)  and  the  intensity  variation  with  fixed  relative  phase.  The  parti- 
cular fixed  relative  phases  A = 0 and  A = ± n are  chosen  here  for 
comparison  because  these  give  extrema  for  Q'  =0,  rr/2,  .t, 

The  effect  of  the  ratio  a2  on  the  extrema  of  first-order  light  intensity 
is  shown  in  Fig.  5-  The  extrema  are  shown  as  a function  of  a2  for  various 
values  of  Q' . Calculations15  of  the  extrema  from  simultaneous  diffraction 
theory  differ  from  those  for  Q'  =0  bv  at  most  the  order  of  10  3 of  the 
total  unit  light  intensity.  However,  the  differences  in  predicted  extrema 
for  larger  Q'  are  shown  here  to  be  significant. 

We  should  remark  that  the  more  interesting  applications  of  successive 
diffraction  theory  are  to  those  cases  for  which  the  results  differ  signifi- 
cantly from  the  case  of  simultaneous  diffraction.  The  simultaneous 
diffraction  equation  equivalent  to  Eq.  (19)  is 

9V  = Z'  Jf  (O  Jq  K)  exp  [-  t (p  dm  + q <JJ] . (25) 

P.i 

Equation  (19)  reduces  to  Eq.  (25)  as  Q and  Q'  approach  zero.  For  suffi- 
ciently small  values  of  %nmpQ  the  argument  of  the  Bessel  function  of 
order  q in  Eq.  (19)  is  negligibly  different  from  v„.  Furthermore,  for 
approximate  equivalence  between  successive  and  simultaneous  diffrac- 
tion, the  argument  of  the  last  exponential  factor  in  Eq.  (19)  must  be 
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small  for  all  values  of  p giving  significant  contributions  to  the  sum. 
We  conclude  that  the  difference  between  successive  and  simultaneous 
diffraction  is  most  marked  for  Q and/or  Q'  large.  But  we  must  recall 
that  complete  validity  of  the  preliminary  Raman-Nath  theories  for 
normal  and  oblique  incidence  has  been  assumed  in  the  development  of 
the  successive  diffraction  theory.  The  limitation  on  the  Raman-Nath 
results  formulated  in  Eq.  (4)  is  closely  related  to  the  parameter  Q by 
the  relationship 

Q^N/v  (26) 

for  the  fundamental  frequency.  This  indicates  that  the  product  Qvn 2 
must  be  less  than  a number  in  the  range  1 to  2,  depending  on  the  accuracy 
required.  An  analogous  limitation  applies  to  the  first  sound  beam. 
In  the  case  of  successive  diffraction,  where  the  final  diffraction  spectrum 
depends  greatly  on  the  amplitude  and  phase  of  the  various  light  compo- 
nents emerging  from  the  first  sound  beam,  the  limitation  expressed  by 

{2nLmAvm)l{n0A*2)^N,  1^iV^2,  (27) 

may  not  be  sufficiently  stringent.  It  remains  to  be  shown  that  there 
exists  a range  of  the  various  parameters  for  which  (1)  the  Raman-Nath 
theories  are  sufficiently  accurate  and  (2)  the  difference  between  successive 
and  simultaneous  diffraction  is  significant.  Reasonable  experimental 
agreement  might  be  obtained  for  Q small  and  arbitrary  Q' , since  Q' 
does  not  affect  validity  of  the  Raman-Nath  theories  but  can  give  signifi- 
cant difference  between  successive  and  simultaneous  diffraction. 
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Introduction. 

Results  of  the  study  of  the  propagation  of  light  in  liquids  under 

pressure  can  be  used  directly  in  the  study  of  the  diffraction  of  light  by 

ultrasonic  waves.  This  fact,  though  it  is  sufficient  justification  for  the 

:tudy,  is  only  one  aspect  of  the  overall  problem,  since  it  might  also  be 

possible  to  gain  basic  information  about  the  structure  of  the  liquid  if  this 

phenomenon  were  completely  understood. 

It  has  been  pointed  out  repeatedly  that  the  equation  given  by 
1 2 

Lorenz  and  Lorentz  , while  It  is  well  founded  theoretically,  is  not  as 

accurate  in  relating  the  index  of  refraction  and  density  of  a liquid  as, 

for  example,  the  EykraanJ  formula,  which  is  a completely  empirical  relation. 

4 

Willard  has  proposed  using  the  Eykraan  formula  to  calculate  the  change  of 
index  of  refraction  with  pressure  in  the  study  of  ultrasonic  waves  in  water 
by  optical  methods.  Although  Willard  did  not  completely  justify  the  use  of 
this  formula,  since  he  referred  to  the  work  of  Gibson  and  Kincaid  who  worked 
only  with  benzene,  it  will  be  seen  that  the  data  collected  here  will  support 
its  use  with  water,  carbon  tetrachloride,  benzene,  and  methyl  alcohol.  These 
data  are  collected  to  show  how  the  values  obtained  by  use  of  the  three 
relationships  agree  with  those  obtained  by  various  experimenters,  and  hence 
to  give  a most  probable  value  of  1/p  dn/dp  to  be  used  in  such  experiments  as 
the  measurement  of  ultrasonic  pressure  amplitude  using  optical  methods. 
Further,  it  is  suggested  that  the  Gladstonc-Dalc  relation  be  used  since  it 
combines  the  advantage  of  simplicity  with  accuracy. 


Discussion  of  Data 

Those  data  are  taken  primarily  from  the  papers  by  Raman  and 

5 

Venkataraman^  by  Gibson  and  Kincaid  , and  Waxier  of  the  Bureau  of  Standards 
who  used  the  Gibson  and  Kincaid  pressure  vessel  and  a Pulfrich  refractometer 
for  making  measurements. 
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In  evaluating  dn/dp,  cognizance  was  taken  of  the  fact  that  the 
curve  relating  n and  p is  not  linear,  and  tnus,  dn/dp  is  a function  of  p, 
in  general  decreasing  with  increasing  p.  Since  the  ultrasonic  measurements 
are  made  at  relatively  low  pressures,  of  the  order  of  one  atmosphere,  an 
attempt  was  made  to  evaluate  the  slope  dn/dp  at  p = 0.  This  was  an  approxi- 
mate procedure  because  of  the  way  the  data  was  taken.  The  data  was  taken  at 
fairly  large  pressure  intervals.  Thus,  it  would  be  expected  that  the  values 
given  would  tend  to  be  smaller  than  the  ones  for  p = 0 because  of  the  curva- 
ture of  the  pressure- index  of  refraction  relationship.  The  typical  curvature 
of  this  relationship  is  shown  in  figure  1,  which  is  a plot  of  some  of  the 
'ata  of  Gibson  and  Kincaid.  From  this  curve  can  be  seen  the  accuracy  of  the 
Jctermination  of  the  slope  at  p = 0.  For  example,  because  of  the  large 
separation  of  the  experimental  points,  the  slope  drawn  on  the  25°  curve  could 
be  1G  percent  lower  if  it  actually  passed  through  the  second  experimental 
point.  One  might  estimate  that  the  accuracy  of  the  slope  is  something 
+ 10  percent.  The  data  for  water  shown  in  Figure  2 from  Waxier  is  more 
linear.  Hence,  the  estimation  of  the  slope  is  more  accurcte.  It  would  be 
still  better  if  the  experimental  points  had  been  taken  at  lower  pressures. 

This  is  the  reason  for  the  fact  that  Raman  and  Venkatareman  used  an  inter- 
ferometric method  for  making  measurements.  Because  of  the  extreme  sensitivity 
of  this  method,  they  were  able  to  make  measurements  when  the  pressure  changes 
were  of  the  order  of  10  centimeters  of  mercury.  They  measured  both  isothermal 
and  adiabatic  values.  Naturally,  the  adiabatic  values  arc  of  more  interest 
if  one  is  studying  the  propagation  of  ultrasonic  waves  using  light  diffraction. 
It  is  to  be  noted,  however,  that  in  general  the  difference  between  the  adia- 
batic and  the  isothermal  values  is  less  than  the  difference  expected  because 
of  experimental  inaccuracy  in  the  ultrasonic  experiments.  For  example,  the 
difference  between  the  adiabatic  and  the  isothermal  values  of  dn/dp  for 
water  as  measured  by  Rarian  and  Venkatareman  is  0.2  percent.  For  carbon 
disulphide  it  is  2 percent,  which  is  the  largest  difference  for  any  liquid 
given  by  them.  In  most  cases,  therefore,  isothermal  values  will  be  sufficient- 
ly accurate  for  ultrasonic  experiments. 

Table  1 is  a compilation  of  the  available  values  of  the  piezo- 
optic  coefficient.  Only  Raman  and  Vonkataranan  give  adiabatic  values.  For 
comparison,  values  are  given  for  both  the  Lorunz-Lorcnfz  equation  and  the 
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Eykman  formula.  Ae  has  been  pointed  out  by  Raman  and  Venice  tar  aman,  the 
Lorenz-Lorentz  equation  gives  values  for  1/0  (dn/dp)  which  are  too  large. 
Without  exception,  the  Eykman  formula  gives  lower  values,  and  almost  without 
exception,  this  lower  value  agrees  better  with  the  experimental  ones.  On 
the  other  hand,  the  Gladstone-Dale  relation  gives  values  which  are  usually 
lower  than  the  experimental  ones  and  which  agree  almost  as  well  as  those 
given  by  the  Eylcman  formulcs.  Therefore,  it  might  be  concluded  that  in  the 
absence  of  experimental  data  the  Gladstone-Dale  relation  is  to  be  preferred. 
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Conclusion 

The  Eykman  formula  gives  a value  for  1/p  (dn/dp)  which  is  of  the 
order  of  6 to  12  percent  lower  than  that  given  by  the  uncorrected  Lorenz- 
Lorentz  equation  for  the  liquids  studied,  and  which  agree  much  better  with  the 
experimental  values.  For  those  liquids  both  the  Eykman  formula  and  the 
Gladstone-Dale  relation  give  values  which  agree  with  the  experimental  values 
of  Raman  and  Venkataraman  to  within  7 percent,  the  Eykman  values  being  in 
general  above,  and  the  Gladstone-Dale  values  in  general  below,  the  experi- 
mental ones.  Listed  in  Table  II  are  the  most  probable  values  for  1/p  (cm/ap). 
The  accuracy  of  these  values  is  certainly  10  percent.  They  are  essentially 
mean  values  of  all  those  found  in  the  literature,  weighted  in  the  direction 
of  the  values  given  by  Raman  and  Venkataraman  since  these  appear  to  be  most 
dependable.  Also  listed  are  the  indices  of  refraction  ond  the  values 
obtained  by  using  the  Gladstone-Dale  relation. 


Table  II. 


most 

Gladstone 

probable 

Dale 

Liquid 

.n 

I fSfi] 

P yapj 

(n-l) 

cci4 

1.458 

0.48 

.458 

Denzcne 

1.498 

0.51 

.498 

Water 

1.322 

0.33 

• 332 

Methyl 

Alcohol 

1.328 

0.33 

.328 
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